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We show how the unitarized description of pion nucleon scattering within 
Heavy Baryon Chiral Perturbation Theory can be considerably improved, by 
a suitable reordering of the expansion over the nucleon mass. Within this 
framework, the A resonance and its associated pole can be recovered from 
the chiral parameters obtained from low-energy determinations. In addition, 
we can obtain a good description of the six S and P wave phase shifts in terms 
of chiral parameters with a natural size and compatible with the Resonance 
^ \ Saturation Hypothesis. 
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I. INTRODUCTION 



The modern way to incorporate chiral symmetry and departures from it in low energy 
hadron dynamics, is by means of Chiral Perturbation Theory (ChPT) [[[[]. Being an effective 
Lagrangian approach all the detailed information on higher energies or underlying micro- 
scopic dynamics is effectively encoded in some low energy coefficients (LEC) which have 
to be determined experimentally. For processes involving only pseudoscalar mesons the ex- 
pansion parameter is p 2 /(47rF) 2 with p their four momentum and F the weak pion decay 
constant [I]. Of course, this expansion works better in the threshold region, breaking down 
at higher energies where the violation of unitarity becomes more and more severe. However, 
it has been shown that this applicability region can be extended by means of unitarization 
methods, describing remarkably well meson-meson scattering and its light resonances up to 
almost 1.2 GeV PH. In addition, it has been shown that the predicted unitarized amplitudes 
can reproduce the threshold region, and provide definite theoretical central values and error 
estimates of the phase-shifts away from threshold and up to about 1 GeV |4]]. 

Baryons can also be included as explicit degrees of freedom if they are treated as heavy 
particles in a covariant framework H called Heavy-Baryon Chiral Perturbation Theory (HBChPT) 
||-||. In this case, the expansion to order iV = 1, 2, 3, ... is written in terms of contributions 
of the form e N /(F 2l M N+1 - 21 ), where I = 1, . . . , [(N + l)/2], and M is the baryon mass. 
The quantity e is a generic parameter with dimensions of energy built up in terms of the 
pseudoscalar momenta and the velocity v M (v 2 = 1 ) and off-shellness k of the baryons. The 

o o 

latter is defined from ps —M v + k, where pb and M are the baryon four momentum and 
its mass at lowest order in HBChPT, respectively. 

Within HBChPT irN scattering, has been calculated up to third order |T0|JTT|] so far[|. In 



general, it is found that the HBChPT convergence is not as good as that of pure ChPT (see 



remarks in Ref. |T(J and also below). Therefore, any attempt to unitarize this amplitude 
based on the increasing smallness of higher order terms, may formally reproduce HBChPT 
series but will likely fail numerically to describe the corresponding phase shifts even in the 
threshold region. Hence, any unitarization method should take into this slow convergence. 

Unitarization methods have also been applied to itN scattering in the literature. Already 
in the seventies Pade approximants had been used to unitarize simple phenomenological 



models [12[], and a more systematic approach based on an effective Lagrangian formalism 
was called for. In addition, several relativistic phenomenological models exist which unitarize 
tree level amplitudes with the K-matrix method providing a reasonably good description of 



txN scattering Hl3[ . More recently, the A and N* resonances have also been considered as 
explicit degrees of freedom within HBChPT (|T1|-|T^]. This requires the introduction of new 
parameters into the Chiral Lagrangian. 

In contrast, the unitarization via the Inverse Amplitude Method (IAM) does not 
introduce new parameters. However, when applied to nN scattering |T7[ , considering the 
higher orders to be increasingly small, the LEC turned out to be very different from those 
found within HBChPT (T(| since they absorb higher order contributions which are not 



negligible. Nevertheless, in [JL7]] it was also shown that the A(1232) can be reproduced with 



1 After submitting this work the HBChPT fourth order result has appeared in the literature [26]. 
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the IAM using 0(p 2 ) parameters constrained to lie within the range of values of Resonance 
Saturation, although the resulting 0(p 3 ) parameters are still very unnatural. This situation 
is somewhat disappointing, since Pade approximants, which are similar to the IAM, together 
with simple phenomenological models, provided good descriptions of ttN scattering. 

In this work we show how the 0(p 3 ) HBChPT series can be reordered and the IAM 
modified, in order to i) implement exact unitarity , ii) comply with HBChPT at threshold 
and iii) describe the A resonance without introducing additional parameters. 



II. PARTIAL WAVE AMPLITUDES 

Customarily, low-energy ttN scattering is described in terms of partial waves ( see, for 
instance Refs P~8| , [TD[ | ) . We will rely heavily on the results in Ref. |TD[, where the third order 



HBChPT expressions where first obtained. However, in that work the full nucleon mass 
dependence has been retained, causing some order mixing. For our purposes, and in order 
to keep track of perturbative unitarity (see below), we have preferred to further expand the 
partial waves in terms of 1/M or 1/F 2 . That is, we only keep the pure 0(e/F 2 ) (first order), 
0(e 2 / (F 2 M)) (second order) as well as the 0(e 3 /F 4 ) and 0(e 3 / (F 2 M 2 )) (third order) terms. 
One thus gets the following expansion, 



which thus coincide with the first, second and third orders in 11 , where the full nucleon 



mass dependence was not retained. However, we also separate these contributions in 

,(i)± _ m ( i,i)±/u;x 
Ii± -pahi [J 

,(2)± _ m 2 (i, 2 )±/^\ 
Ji± - F 2 M h± [J 

where tf^ 71 are dimensionless functions of the dimensionless variable u/m, independent of 
F, M and m. The analytical expressions are too long to be displayed here, but can be easily 
obtained starting from Ref. . The convenience of the double superscript notation will be 
explained below; n + 1 indicates the power in 1/F and m the total order in the HBChPT 
counting. The unitarity condition Im/^ 1 = — q becomes in perturbation theory 

Im f P/> ± = Im t?f ± = Im ± = (3) 
Imt^^h 2 (4) 

The equivalent of the last equation above for the / amplitudes is not satisfied exactly if 
the nucleon mass dependence is retained through a redefinition of the nucleon field as in 



Ref.( |T0|) (see, for instance, comments in [] fj and [|17[), although, of course, the corrections 



are just higher order in HBChPT. These exact relations for the perturbative contributions 
will be very convenient later on, and, as we have commented above, that is the reason why 
we have preferred to expand the amplitudes as in Eq. ([I]). 
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The scattering lengths a^± and effective ranges bf ± are defined by 

Ref±=q 2l (a±+q 2 b±+-- 



(5) 



, a l ± and b l ± , 



Obviously, the expansion of Eq. (g) carries over to the threshold parameters 
which we will use next to illustrate the slow low energy convergence of the HBChPT series. 

Throughout this paper we use F = 92.4 MeV, M = 938.27 MeV, m = 139.57 MeV and 
qa = 1-26. Concerning the LEC , there are several determinations: The first one, which we 
give as Set I in Table I, was obtained from a fit to the extrapolated threshold parameters a^ + , 
a t±' the Goldberger-Treiman discrepancy and the nucleon a— term |IJ. Note that 



the full nucleon mass dependence was retained. For our purposes it is thus more convenient 
a more recent determination |TT| from a low energy fit to ttN phase-shift data, although the 
authors used a different notation for the chiral parameters. The resulting LEC translated to 

19 11]]) are given in the Set II of Table I. On the theoretical side, 



the cij and hi notation (see 
there are estimations of the a, parameters assuming they are saturated by the exchange of 
resonances ||, in fairly good agreement with experimental determinations. 



III. LOW ENERGY CONVERGENCE 

In all our following considerations we further expand the amplitudes and threshold pa- 
rameters as in Eqs.(@). For the scattering lengths we thus have 

4, 2 X M = + ^ M + + f> 3 > + • • • (6) 

and a similar expression for b\j 2 jm 2l+3 . The slower convergence of the HBChPT series 
as compared with ChPT was pointed out by the authors of the first calculations, (see for 
instance the comments in |TT[] ) . For our purposes, we have found instructive to separate 
the contributions to the scattering lengths and effective ranges of the lowest partial waves, 
which are given in Table II, using the set II of LEC in Table I. Note that, in this way, the 
threshold parameters are predictions, unlike those from set I which uses them as an input. 
A distinctive pattern emerging from this Table is that the contribution of order 1/(F 2 M 2 ) 
is always rather small. Only in some cases, however, is the contribution of order 1/F 4 also 
small. This is so in the P33 channel in particular, which is the one less likely to be well 
described by HBChPT alone due to the low mass of the A(1232) resonance. Hence, it 
seems that close to threshold the 1/F 2 expansion converges faster than the 1/M expansion. 
Actually, the 1/F 2 and 1 / (F 2 M) terms are comparable. It is clear that any unitarization 
method will only be consistent with the HBChPT approach, if it treats both the first and 
the second order as equally important. 



IV. UNITARIZATION METHOD FOR THE REORDERED SERIES 

Our unitarization method assumes that, as suggested by the perturbative calculation, the 
chiral expansion in terms of 1/F 2 converges much faster than the finite nucleon mass 1/M 
corrections. Indeed there are some recent theoretical attempts [22,23] to define a relativistic 
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power counting not requiring the heavy baryon idea, returning somehow to the spirit of older 



relativistic studies f24| . As a matter of fact, we show that the well-known IAM approach 
applied to the reordered HBChPT series generates the P 33 phase-shift satisfactorily, just 
using the LEC determined at low energies. In addition, it is possible to improve the overall 
description of the six ttN scattering S and P waves, using those very same parameters. 
Finally, if one fits the ttN phase shifts with this method (either constraining the parameters 
to the Resonance Saturation Hypothesis or leaving all the parameters unconstrained) the 
resulting parameters have a much more natural size, and the x 2 P er d.o.f is considerably 
better than those obtained with the IAM applied to plain HBChPT |T7[ ]. 





M.Mojzis @ 
Set I 


N.Fettes et al Q 
Set II 


Resonance Saturation fit 
Set III 


Unconstrained fit 
Set IV 


dl 


- 2.60 ± 0.03 


- 2.69 ± 0.4 


-2.065 ± 0.007 


- 1.36 ± 0.02 




1.40 ± 0.05 


1.34±0.1 


0.915±0.005 


0.438±0.015 


a3 


-1.00 ± 0.06 


-1.15±0.1 


-0.85 (input) 


-0.70±0.04 


a 5 


3.30±0.05 


3.1±0.5 


2.700±0.001 


1.29±0.04 


61+62 


2.40±0.3 


2.60±0.2 


3.95±0.04 


3.06±0.3 


63 


-2.8±0.6 


-3.96±0.9 


-1.45±0.03 


-0.41±0.27 




1.4±0.3 


0.55±0.5 


-1.17±0.17 


-1.5±0.2 


bw ~ h5 


6.1±0.6 


7.1 ±0.5 


5.86 ±0.19 


7.4 ±0.5 


big 


-2.4±0.4 


-1.72±0.3 


-0.44±0.21 


-3.7±0.2 



TABLE I. HBChPT low energy constants. Those in the first column were obtained from fitting the 
extrapolated threshold parameters a^ + , b$ + , af_j_, the Goldberger-Treiman discrepancy and the nucleon 
a-— term to the HBChPT predictions ]io| . In the second column we give the parameters obtained from a low 
energy fit to ttN phase-shift data pT| . In the third column we give the parameters obtained from an IAM 
fit with the ai constrained to the ranges predicted by resonance saturation. Finally, the fourth column is 
the result of a totally unconstrained fit. 





l/F 2 


1/(F 2 M) 


1/(F 2 M 2 ) 


l/F 4 


Total 


SP98 


KA85 




-0.65 


-0.04 


±0.07 


±0.07 


-0 ^+ u - ie 


-0.64 ± 0.01 


-0.72 


<i 


1.3 


-0.34 


-0.08 


0.07 


0.94±0.23 


1.27 ± 0.02 


1.26 


a 3 3 


35.3 


47.95 


-1.75 


0.26 


81.818J 


80.3±0.6 


78.75 


ai 3 


-17.7 


15.46 


-3.1 


-6.35 


-11.66±0.9 


-10.5 ± 0.9 


-11.00 


031 


-17.7 


12.96 


-1.84 


-9.77 


-16.332 


-15.9 ± 1.0 


-16.19 




-70.7 


85.6 


-3.77 


-46.13 


-35.0l!J 


-27.0 ±1.4 


-28.67 



TABLE II. HBChPT, lowest partial S and P wave- scattering lengths, a l 2I2J (in GeV _2i_1 units) for 
ttN scattering using set II in Table I, decomposed as a sum of terms of first order l/F 2 , second order, 
1/(F 2 M), and third order, 1/(F 2 M 2 ) and 1/(F 4 ). The sum of all the terms yields the total scattering 
length parameter. For brevity, we only quote the errors in the final sum. The experimental values come 
from Ref. |21| for SP98, and Ref. [^0| for KA85. For very recent and accurate values of only the S-wave 
scattering lengths see Effl 
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The IAM is based on the fact that elastic unitarity imposes the following relation on a 
generic partial wave / and its inverse / _1 (we drop the I, I and J labels for simplicity): 



Im/ = g|/| 2 => Imf-^-q (7) 
As a consequence, any amplitude satisfying exactly elastic unitarity has the form 

f = ReA-zq (8) 

Thus, we only have to calculate Re/ -1 , whose different approximations provide different 
unitarization methods. Here we consider its expansion in terms of m 2 /F 2 , i.e. 

3 

f(u,m,F,M) = ^t {1 \uj/m, m/M) + ^-t {3 \uj/m, m/M) + ... (9) 
F 2 F A 

where we have used that p 2n+1 ' are dimensionless functions, only depending on dimensionless 
variables. The functions and are only known in a further m/M expansion, 

t^ 2n+l \uj/m,m/M) = t^ 2n+1 ' 2n+1 \uj/m) + H t {2n+1 > 2n+2 \u / m) + Q^V t^ 2n+l ' 2n+3 \u /m) . 

(10) 

yielding Eq. (P and Eq. (0) after a suitable isospin projection^. 
Perturbative unitarity in this expansion requires, 

lmt« = , lmt® = ±\tW\ 2 , (11) 

m 

which also imply an infinite conditions in the 1/M expansion. For / _1 we get then 

1 F 2 1 



m- 



f m tW [tW] 



+ ... (12) 



Obviously, expanding in m/F 2 may be justified provided these corrections are small. As 
we have shown, for the P33 channel, they are small precisely at threshold, and there the 
unitarization scheme will, approximately, reproduce the perturbative result. At the same 
time, unitarity is exactly implemented since, thanks to eqs.flTTD the above formula is of 
the form given in eq.(H). However, the 1/M terms are not small corrections at threshold, 
and thus we do not further "expand the denominator" . In this way we keep the first mass 
corrections as equally important. At present, only the t^ 1 ' 1 ', t^ l,2 \ t^' 3 ' and t^ 3 ' 3 ^ HBChPT 
terms are known (see our previous footnote [I]). Therefore, we do not know neither t^ 3 ' 4 ^ 
nor f( 3 ' 5 \ i.e 0(1/(F 4 M)) and 0(1/(F 4 M 2 )) respectively, and we can only approximate the 
numerator of the second term by ~ t^ 3 ' 3 \ That is, we know the first three orders of 
the denominator l/m expansion, but not their counterparts in the numerator. Therefore, 



The relation between f l 2I2J and ff ± is given by fi 2l±1 = f+± - f l± and fi 2l±1 = f+± + 2f l ± . 
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although we could use t^ 1 ) ~ t^ 1 - 1 ) + t^' 2 ' +t^ 1 ' 3 ^ in the denominator, for consistency with the 
numerator expansion we keep only t' 1 ' ~ t^ 1 ' 1 '. Incidentally, this corresponds to the static 
limit of the second term. Note also that strict unitarity is still satisfied. Of course, it would 
be desirable to compute at least t^ 3 ' 4 ^ and t^ 3 '^ in order to be able to keep also t^. 2 ) in the 
denominator of the second term of the mentioned equation. After these remarks, we have 

1 F 2 1 

m 7777^ ( 13 ) 



/[unitarized ™ + § t^ 2 ) + (f ) 2 [t™] 

At threshold, our formula yields a modified scattering length 

1 F 2 1 c/ 3 ' 3 ) 

" ™7T7TTm ( 14 ) 



aunitarized Tfl C^ 1 ' 1 ) + % ftW) + (§) 2 ad' 3 ) [c^' 1 )] 2 



1 a-C 3 - 3 ) 

m T77Tn^ ( 15 ) 



a HB chPT - (f) «( 3 ' 3 ) ""[ad- 1 )] 

which, using set II in Table I, yields CI33I unitarized = 82.95^o;35GeV~ 3 , to be compared with 
dgglnBChPT — 81.8lQ;gGeV -3 , both compatible with the experimental values. Notice that if 
we had expanded the denominator considering the second order contribution 0(1/(F 2 M)) 
to be small we would have obtained strictly the IAM method as used in Ref. p 



OlAM 




m a (1 > 2 ) /m\ 2 a (1 > 3 ) /m\ 2 [a (1 > 2 )] 2 l a^ 3 ' 3 ) 

m [«(y)] 2 ~ vm/ [ a (i,D]2 + (.mJ j a (i,D]3 j ~ m [ a (i,D]2 ^ 16 ^ 



yielding ci^Iiam = 22.8GeV 3 . This explains why an unconstrained IAM fit with standard 
HBChPT leads to LEC which are so different from those found in [10], as noted in Ref. [17 



V. NUMERICAL RESULTS 

From our previous discussion it is clear that at threshold our unitarized amplitude will 
reproduce very accurately and within error bars the HBChPT results and hence the exper- 
imental data. In addition, we expect that the P33 phase shift can be extended up to the 
resonance region by propagating the errors of the LEC. It is also tempting to extend the 
other five S and P wave phase shifts. We show in Fig.l the phase shifts obtained from our 



unitarization method, Eq. (0), compared with the experimental tcN data pl| . The shaded 
area corresponds to the phase shifts obtained by propagating the errors of the parameters 
given in Ref. |TT[| (See Table I, set II) by means of a Monte Carlo gaussian sampling of the 



LEC for any given CM energy value. Only for comparison, the dotted line corresponds to 
standard HBChPT extrapolated to high energies. 

As one can see from the figures, the prediction of our unitarized approach produces a 
distinctive resonance in the P33 channel, with very similar parameters to the physical A as 
we will see below. Concerning the other channels, there is some improvement in the S waves, 
and a worse behavior for the P13, P31 and Pq, but note that these three partial waves have 
very tiny phase shifts, and any small error yields a large relative deviation. 

The mass and width of the A resonance can be obtained either from the phase shifts, by 
means of 5l 3 \ s=M ^ = n/2 and 1/Ta = M&(ddl 3 /ds)\ s=M 2^, or from its associated pole in the 
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FIG. 1. Phase shifts as a function of the total CM energy ^/s. Experimental data are from Ref. 
The shaded areas correspond to the propagated errors of the parameters in set II, which were obtained from 
low-energy data. The dotted line is the HBChPT result extrapolated to high energies. The dashed line 
is a fit with our unitarization procedure constrained to Resonance Saturation, whereas the continuous line 
corresponds to an unconstrained fit. 

second Riemann sheet (y/s po i e ~ Ma — iT^/2). We give in Table III the results for different 
parameter sets. Note that the width of the "predicted" resonance from the parameters 
determined from low-energy data (set II) is qualitatively very similar to the real A. Of 
course, once we fit to the data (set III and IV), we obtain a much better description. 

However, it was pointed out in Ref. jl7| that the direct fit using the IAM directly on 
the HBChPT series leads to chiral parameters of very unnatural size. That is not the case 
when we fit with the reordered method proposed here as it can be seen in set IV of table I. 
This set comes from an unconstrained fit to the six S and P wave irN phase shifts, which is 
represented as a continuous line in Fig.l. For the fits, which start at 1130 MeV, we have used 
the MINUIT minimization routine assigning a 3% uncertainty as in [|TIJ] plus a systematic 
error of one degree to the data in |H] (similar treatments are followed in |Ilj|{P7j| ). 
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Much more interesting are those fits where the 0(p 2 ) parameters are constrained to 
the range estimated by the Resonance Saturation Hypothesis || . The fitted parameters are 
given as set III in Table I, and the result is represented as the dashed line in Fig.l. In this 
case the P n and Sn are not so well described, which may be due to effects of the iV*(i440) 
and iV(1535), respectively, which are the closest resonances to the energy regions displayed 
in Fig.l. Indeed, the former plays a marginal role in the Resonance Saturation Hypothesis 
whereas the latter is not even considered. 

A particularly relevant feature of these fits is that not only the resulting parameters have 
a more natural size, but also the x 2 P er d.o.f. is between three and four times smaller than 
for the corresponding IAM fit applied to the standard HBChPT ordering. From this we can 
conclude that considering the 1/M expansion separately as in our formalism is a sensible 
approach, apart from the details of its precise realization. 
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FIG. 2. P33 phase shift. The upper shaded area corresponds to the result of propagating the errors of the 
set I parameters. This illustrates the uncertainties due to the choice of parameter sets from the literature. 
The intermediate shaded area cover the propagated uncertainties of set II if we applied the IAM in the static 
limit. Finally, with the same set we show, in the lowest shaded area the result of the IAM when applied 
directly to the standard HBChPT expansion. 





set I 


set II 


set III 


set IV 


PDG 1 25] 


M A (MeV) 


1240 tf 7 


1222 t\l 


1227 ±4 


1226+{^ 


1230 - 1234 


T A (MeV) 




117+24 

11 '-18 


104.1 til 


W7tH 


115-125 


Re (pole position) (MeV) 


1205 


1204 


1204 


1204 


1209-1211 


T -2 Im (pole position) (MeV) 


110 


110 


110 


84 


98-102 



TABLE III. A resonance parameters. The two first rows are obtained from the condition 
^33 1 s=M^ = 77 /2 and 1/Ta = Mj\{d5\^/ ds)\ s=M 2 . The second two rows are the real part and 
minus twice the imaginary part of its associated pole position, calculated for the central values of 
each set (the errors are expected to be of the same order than for the first two columns). 



S 



As an illustration of the uncertainties due to the different determinations of chiral param- 
eters, we have also shown in Fig. 2 the area corresponding to our formula applied to set I. We 
reobtain the same qualitative result, although numerically the mass and width of the A are 
worse than those obtained with set II. In addition, in order to estimate the convergence rate 
of our calculation, we have also plotted in Fig.2 the prediction in the static limit (M — > oo). 
The shaded area in Fig.2 corresponds to the propagated errors of the parameter set II in 
this limit. As we see, there is also a distinctive resonant behavior, so that the bulk of the 
dynamics is contained in the static limit. However, the finite mass corrections, particularly 
the 1/(F 2 M) contribution, are important to achieve a better description?]. 

We have also studied what happens if one includes the incomplete higher order contribu- 
tions to the second quotient in Eq. (0), i.e. if one approximates ~ t^ 1 ' 1 ' + ^ 1,2 ^ + ^ 1,3 ) in 
its denominator. In such case we obtain a worse result, closer to that of the static limit, but 
still there is a distinctive resonant behavior, improving the IAM results with the standard 
HBChPT expansion. This suggests that the unknown higher order contributions (see foot- 
note |l|) in the numerator could give rise to some cancellation with those still incomplete of 
the denominator. Finally, we also show in Fig.2 the results obtained within the conventional 
IAM approach, for the parameter set II. As it was already pointed out in |17| the result is 
extremely poor if one uses the parameters determined from low energy data. 



VI. CONCLUSIONS AND OUTLOOK 



Heavy Baryon Chiral Perturbation Theory provides definite predictions for the irN scat- 
tering amplitudes in the threshold region. However it violates exact unitarity if the per- 
turbative expansion is truncated to some finite order and also is unable to describe the A 
resonance (and its associated pole) in the P33 channel. The analysis up to third order shows 
that the leading finite nucleon mass correction, which is second order, is of comparable size 
to the static approximation and in fact it dominates the corrections at threshold. This sug- 
gests a unitarization method using the expansion in inverse powers of the weak pion decay 
constant but without making the heavy baryon expansion. Such an idea is supported by 
recent theoretical attempts to redefine a relativistic chiral counting for baryons. We have 
proposed a unitarization scheme based on the Inverse Amplitude Method applied to this 
reordered HBChPT expansion. It provides a prediction for the irN phase shifts, which gen- 
erates a A resonance from the low energy constants and their errors, as determined from 
HBChPT. The fits within this scheme provide chiral parameters of a natural size and a bet- 
ter overall description than those performed with the IAM applied to the HBChPT standard 
expansion. This result suggests that including the 1/M expansion separately is a sensible 
physical approach. In addition, this method can be easily generalized to higher orders and 
coupled channels. Further work along these lines is in progress. 



3 The 1/(F 2 M 2 ) correction turns out to be quite small. That was expected, since it neither 
provides a sizeable contribution at threshold (unlike the 1/(F 2 M) correction ) nor it is responsible 
for the restoration of unitarity (like the 1/F A correction). 
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